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It has been suggested that the supermassive black holes, at the centers 
of galaxies and quasars, may initially form in single collapses of relativistic 
star clusters or supermassive stars built-up during the evolution of dense star 
clusters. We show that it may be possible for ICECUBE (a planned 1 km^ 
neutrino detector in Antarctica) to detect the neutrino bursts associated with 
those collapses at redshift z < 0.2 with a rate of 0.1 to 1 burst per year. 
Such detections could give new insights into the formation of structure in 
the universe, especially when correlated with gravitational wave signatures or 
even gamma-ray bursts. 
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In this letter we describe a possible means for detecting the neutrino bursts accompanying 
the formation of supermassive black holes at cosmological distances. (Hereafter the term 
"supermassive" implies a mass > 3x 10'' Mq.) Mounting evidence suggests that supermassive 
black holes are fairly common in the universe. Observations with the Hubble Space Telescope 
point to supermassive black holes in nearly every galaxy examined so far Likewise 
supermassive black holes have long been thought to be the central engines of quasars and 
Active Galactic Nuclei (AGN) 0. The masses inferred for these black holes typically range 
from ~ 10^ Mq (as in our galaxy) to ~ 10^ Mq (as in some quasars). It is conceivable that 
every galaxy, large or small, harbors a supermassive black hole at its center. In other words, 
supermassive black holes may have an abundance of one per IO^^Mq (baryon mass) object, 
the typical baryon content of dwarf galaxies such as M32. Such a conjecture is entirely 
consistent with observations if, like those nearby ones, most galactic nuclei are inactive P]. 

How these supermassive black holes were formed is still a mystery. However, a natural 
formation route is through accretion onto a seed black hole. The seed black hole presumably 
must be formed earlier from the collapse of a single object or a cluster of objects. It is likely 
that these seed black holes were themselves supermassive, having masses perhaps only an 
order of magnitude or two below that of the current end-state black holes. This is because 
it would be very difficult for a much smaller seed black hole, such as that resulting from 
an ordinary solar-mass-scale stellar collapse, to accrete material efficiently enough so as to 
reach a supermassive mass scale in a Hubble time at an early epoch @,||. 

It has been suggested that the supermassive seed black holes may form from collapses of 
dense star clusters or supermassive stars that were built up as an intermediate phase during 
the final evolutionary stage of collapsing star clusters 0. Dense star clusters have been 
observed to reside at the center of galaxies and in general have been invoked to explain 
AGN activities The formation of an intermediate-stage supermassive star has been 

argued to be on several of the possible evolution routes of collapsing star clusters 0. 

A supermassive star with mass M > 3 x 10*^ Mq is an index n = 3 polytropic config- 
uration with a high entropy per baryon S ~ 300 (M/IO^Mq)^/^ (in units of Boltzmann's 
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constant) [T^. A supermassive star eventually collapses into a black hole as a result of 
the Feynman-Chandrasekhar |Tl| instability. Numerical calculations also have shown that a 
dense star cluster of compact objects would eventually collapse as a result of the Feynman- 
Chandrasekhar instability just as would a single supermassive object [|12|. 

If a supermassive black hole formed as a result of a single collapse, whether from a single 
supermassive star, or a supermassive relativistic cluster of compact objects, we would expect 
the collapse to be accompanied by a burst of thermally-produced neutrinos. The neutrino 
emission would carry away a significant fraction of the gravitational binding energy of the 
collapse, which is of order the rest mass energy of the homologous core, ~ erg 
(where homologous core mass, denoted by {10^Mq)M^^, is the mass that plunges through 
the event horizon as a unit). The homologous core mass M|^^ is only a fraction of the initial 
stellar mass M™'* = M^'^/IO^M©, with the fraction determined by the entropy loss during 
the collapse |T^. This fraction is roughly ~10% for a non-rotational and non-magnetized 



spherical supermassive star of M™'* ~ 10. 

If the homologous core is massive enough (namely M^'^ > 0.1), it will be transparent to 



neutrinos [|T^. This is in stark contrast to the case of an ordinary core collapse supernova, 
where the neutrinos are trapped in the core. Since in our case the collapsing material will 
be essentially in free fall, we should have a set of well defined neutrino luminosity /time tem- 
plates with which to compare to observations. Such time templates of neutrino luminosity 
have a characteristic gradual rise followed by a sudden drop due to the core becoming a black 
hole. Figure 1 shows an example of such a time profile for the (un-normalized) neutrino lu- 
minosity, and the time profile of the corresponding (un-normalized) signal in ICECUBE. 
The duration of the neutrino burst is typically < tdyn, where tdyn is the dynamic collapse 
timescale. For a non-rotating, non- magnetized progenitor, tdyn ~ M^'^ seconds. 

Shi and Fuller |T3| showed that for the collapse of a non-rotating non-magnetized spher- 
ical supermassive object the neutrino burst could carry away a fraction ~ 0.04(M|^*")~^'^ 
of the total gravitational binding energy (assuming M^'^ > 0.1, so that this fraction is less 
than 1). The neutrino release is equally partitioned between neutrinos and anti-neutrinos. 



with v^Ve accounting for 70% of the total neutrino flux. The average neutrino energy is 
similar for all species and is ~ 4(M|^*-^)~°'^ MeV. 

If the progenitor is rapidly rotating or if magnetic stresses are appreciable, the collapse 
timescale could increase significantly. Neutrinos would then have a much better chance of 
escaping before the core moves through an event horizon. In such a case, the neutrino fiuence 
could be an order of magnitude higher than in the non-rotating case (but still limited by the 
total gravitational binding energy), and the average neutrino energy could be a factor of 2 



higher as well ||13[ . The partition of energy among the different neutrino species, however, 
would remain the same. 

Though the neutrino emission accompanying the collapse of a supermassive object is 



gigantic, it is almost impossible to detect if it originates at a redshift z > \ ||T3|. But 



such a neutrino burst may become detectable if the collapse occurs at a redshift z < 0.2. 
These recent supermassive collapse events may not be uncommon, since there are still large 
numbers of quasars (and AGNs) at redshift < 0.2. 

To estimate how many supermassive progenitors might collapse at a redshift z < 0.2, we 
need to know how abundant they are and how their collapse rate evolves with redshift. A 
natural and reasonable starting point is to assume an abundance of one per galaxy (including 
dwarf galaxies, i.e., one per ~ lO^^M© of baryons), and assume a collapse rate that traces 
the formation rate of quasars, because quasars are powered by supermassive black holes. 
If we assume that the average quasar lifetime, which is short (~ 10^ years) compared to 
the Hubble time, does not evolve with redshift, the quasar formation rate would then trace 
quasar number density as a function of redshift. In figure 2 we replot the evolution of the 
quasar number density (proper) as measured by Shaver et al. [|15|, as well as that inferred 
from the quasar catalogue compiled by Veron-Cetty and Veron |T^. The plot assumes 



a Hubble constant Hq = 50 km/sec/Mpc and a deceleration parameter go = 0.5, which 
will be the working assumption in the following calculations. (The details of cosmology do 
not affect our results significantly.) The sample of Shaver et al. is radio- selected and is 
argued to have the least amount of selection bias ||l5l. We thus adopt the quasar number 



density evolution in the redshift range z > 0.5 as implied by this sample. The catalogue 
of Veron-Cetty and Veron, on the other hand, is a compilation of quasars brighter than 



absolute magnitude Mb = —23 from different samples with different selection criterion |]T6 
Its degree of completeness is uncertain. Therefore, even though the latter catalogue has 
far more data at the low redshift end, its inferred quasar number density evolution below 
z < 0.3 cannot simply be taken as the extension of the Shaver et ai.'s evolution function 

z < 0.3. The catalogue does, however, seem to give a number density evolution in fair 
agreement with Shaver et al.^s evolution function at 0.5 < z < 2. If nothing else, the 
Veron-Cetty and Veron catalogue may indicate that the quasar number density at z < 0.3 
potentially is comparable to the quasar number density at 2; ~ 0.5. To account for the 
uncertainty in quasar number density evolution at the very low redshift, we will calculate 
the nearby supermassive black hole formation rate for three values of the relative quasar 
number density in the redshift range 0.1 < z < 0.2 (the three thick solid lines in figure 2), 
corresponding to a quasar number density that is the same as (case 1), one-third of (case 
2), and one-tenth of (case 3) the quasar number density at 2; ~ 0.5. 

We denote the collapse rate of supermassive objects as 0, whose redshift evolution is 
inferred from figure 2. Assuming a matter-dominated flat universe, we find that the total 
baryon mass in a unit volume that was incorporated into supermassive black holes is now 

/■*o (A M 

P-^= n \ ,3 dt^ 0.02500 M to, (1) 

where to is the age of the universe today, M is the mass of the supermassive black hole 
today, and where 0o is the normalization required to give the true formation rate in a 
unit volume. On the other hand, we have assumed psm = {M / 10^^ Mq) pi, where pb ~ 
1.88 X 10~^^fib^^ g cm~^ (with Qb being the baryon density in units of the critical density 
and h = ifo/(100km/sec/Mpc) is the baryon density today. Therefore, 

,„...10-"(I^i|)(l^)Mpc-.- (2, 

The event rate as observed now of these supermassive collapses that occurred in the 
redshift range zi < z < Z2 is 



where r is the comoving spatial coordinate. Note that in this equation, in addition to the 
factor (1 + z)~^ stemming from volume expansion, there is an additional factor (1 + z)~^ 
due to the cosmic time dilation. In a matter-dominated universe with go = 0.5, we have 



r = 3cto(l — l/vT+^), where c is the speed of light. The rate of those collapses at 
0.1 < 2 < 0.2 is then 

where Rq ~ 1, 0.3, and 0.1, for case 1, 2 and 3, respectively. (The rate for all supermassive 
collapses is roughtly 0.1 per day with the same scaling factors.) These rates scale inversely 
with the amount of baryon mass that contains one supermassive black hole, which we have 
assumed to be IO^^Mq on average. While we are certainly very interested in collapses at 
z < 0.1, we cut off the collapse rate at z = 0.1 to avoid extrapolating the quasar number 
density evolution function too far toward low z. Even so, the rate of nearby collapse of 
supermassive black hole progenitors is potentially much higher than the Type II (and Type 
lb) supernova rate in our galaxy. 

With these estimated characteristics, the neutrino bursts from such relatively nearby 
supermassive object collapse events may be detectable with ICECUBE, a ~ 1 km^ neutrino 
detector with ~ 10^ optical modules soon to be proposed. Just like the currently-operating 
AMANDA (Antarctic Muon and Neutrino Detector Array) detector, it will have sensitivity 
to bursts of neutrinos with energies > 10 MeV [jl^. The principle is similar to that of 



detecting supernova neutrinos, i.e., to detecting the Cerenkov photon flashes resulting from 
relativistic positrons produced by the reaction + p n + . The cross section for this 
process is ~ 9 x 10~^^(i?p^/l MeV)^ cm^. The Cerenkov photons are collected by optical 
modules and their numbers are counted. Using the module efficiency and threshold estimates 

we 



in Halzen, Jacobsen and Zas [jT^ , and the Ve spectrum calculated by Shi and Fuller |]T3 
find that the number of events expected in ICECUBE from a neutrino burst resulting from 
the collapse of a supermassive object would be 



~ 10 = iV„ (— ^) ( ,,i,.,(^HC)-o..,,g ) « (5) 

where Nm is the number of optical modules in the detector, d is the proper distance to the 
source of the burst {d ^ 750 Mpc for z = 0.15), E^p is the total energy release in neutrino 
emission, and a ~ C(l) accounts for the energy spectrum of the Pg emission. The expected 
number of events per optical module (iVevent/-^M) from collapses of homologous cores 
with various masses ai z = 0.15 is shown in Table 1. In our calculation, we have imposed 
a 30% artificial limit on the fraction of the total gravitational binding energy that can be 
carried away by neutrino emission. We believe this fraction is attainable, considering that 
much higher efficiency is possible by magnetic mechanisms. 

Optimal neutrino signal outputs are obtained from collapses with M^*-' ~ C(O.l) 
and with rotation/magnetic fields. For these conditions the average neutrino energy 
(cx (M|^*-')~°-^) is high, as is the fraction of gravitational binding energy carried away by neu- 
trinos. Rotation and/or magnetic fields are always likely to be present during the collapses. 
It is also possible that the lower mass supermassive progenitors that yield M|^'-' = 0(0.1) 
dominate the population of supermassive objects. The optimal conditions of neutrino de- 
tection, therefore, may in fact represent most of the supermassive collapses. 

Competing with the neutrino signal is the fiuctuation of the background level of ICE- 
CUBE, which at the current AMANDA level should have a standard deviation of 20A^^^ 
counts in 1 second. To obtain a signal-to-noise ratio (S/N) of unity for a collapse at 2; = 0.15 
with M|^'-' = 0.3 would then require 2.8 x 10^ optical modules. If we take advantage of the 
time template of the neutrino signal, we may not require a signal to noise as high as S/N= 1. 
A S/N of 0.1, for example, would only require 2.8 x 10^ optical modules for a detection. The 
rate of such detections could be of the order of 0.1 to 1 per year, if most of the supermassive 
collapses at z < 0.2 have rotating/magnetized homologous cores with relatively lower masses 
in the supermassive mass spectrum, M^'^ — 0(0. 1). 

Since to maintain the same S/N, the required Nm is proportional to d^, a chance detection 
of a supermassive collapse at z < 0.1 (and hence d < 500 Mpc) would make the detection 
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scheme even more attractive. Note also that the number of modules (Nm) required scales 
linearly with the duration of the neutrino burst, which we have taken to be tdyn in the 
above estimates. If this duration is only a fraction of tdyn as shown in the example in figure 
1, the number of optical modules can be further reduced. In addition, if the energy of 
neutrino emission during the collapses can exceed 30% of the total gravitational binding 
energy, Nm will scale down as the inverse square of this fraction. Therefore, the detection 
potential actually starts at the Nm ~ 5 x 10^ level, which is currently envisaged for the 
expansion of the operating AMANDA neutrino telescope to the ICECUBE configuration. 
An ICECUBE detector with the order of 10^ optical modules has already been contemplated 
for other science missions, e.g., for searching for nucleon decay with a sensitivity that cannot 
be matched by conventional techniques. 

The above calculations have been based on the detection technology available to 
AMANDA, without any change in design of the planned ICECUBE. This leaves room for 
future improvement. Moreover, if the formation of supermassive black holes via supermas- 
sive object collapse also gives rise to 7-ray bursts W^, or gravitational radiation detectable 
in proposed low frequency gravitational wave detectors such as LISA, we could enhance 
our ability to detect the associated neutrino bursts by knowing when to look. In addition, 
the combined neutrino/7-ray burst /gravitational wave signal would offer a golden chance to 
explore many questions in particle physics, astrophysics and gravitational physics. 
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TABLE I. Expected Neutrino Events per Optical Module in ICECUBE from Collapses at 

z = 0.15. 





Rotation/Magnetic Field 


OL 




TVc-irc-n-f / TV A/f 


0.3 


No 


0.8 


7% 


8 X 10-^ 


0.3 


Yes 


9 


30%^ 


1.2 X 10-2 


1 


No 


0.01 


4% 


1 X 10"^ 


1 


Yes 


1.2 


30%" 


5 X 10"^ 


3 


Yes 


0.017 


22% 


1.7 X 10-^ 



^Saturating our imposed limit. 
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Figure Captions: 

Figure 1. Solid line: the time profile of the (un-normalized) neutrino luminosity from the 
collapse of a non-rotating non-magnetized spherical supermassive object [|13[. Dashed line: 
the time profile of the neutrino signal in the ICECUBE detector. The time axis is in units 
of tdyn, the dynamic collapse timescale. 

Figure 2. The relative proper number density of quasars as a function of redshift, normalized 
to peak at 1. Data points are the measurements of Shaver et al. [^]. The thin solid line 



is the fit of Shaver et al. ITq] to the data. The dashed line is inferred from the quasar 



catalogue of Veron-Cetty and Veron ||I6[, without any correction for bias. The three thick 
lines represent our adopted values for the relative quasar number density in the redshift 
range 0.1 < z < 0.2. 
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